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ABSTRACT 

During the tidal disruption of a main sequence star by a massive black hole (BH) hav- 
ing mass Mbh ^ 10 7 M Q , the stellar debris is expected to fall back to the BH at a rate 
well above the Eddington rate. Some fraction of this gas is predicted to be blown away 
from the BH, producing an optically bright flare of radiation. We predict the spectra 
and spectral evolution of tidal disruption events, focusing on the signatures produced 
by photoionized gas outside the photosphere of this super-Eddington outflow. We show 
that the spectrum of such an outflow should show absorption lines that are strongly 
blueshifted relative to the host galaxy, are typically very broad (0.01 — 0.1c), and are 
most prominent at ultraviolet wavelengths (e.g., C IV, Lyman a, O VI) at early times 
(<; 1 month for a ~ 10 6 Mq BH). There may also be optical absorption lines of hy- 
drogen and He II if there is a lower velocity component to the outflow (<; 0.01 c). At 
later times, the outflow falls out of thermal equilibrium and the continuum emission 
likely hardens — the absorption lines will then primarily be in the soft X-rays. 

Supernovae in galactic nuclei are a significant source of confusion in optical surveys 
for tidal disruption events: we estimate that nuclear Type la supernovae are two orders 
of magnitude more common than tidal disruption events at z ~ 0.1 for ground-based 
surveys. Nuclear Type II supernovae occur at a comparable rate but can be excluded by 
pre-selecting red galaxies. The contamination from nuclear supernovae can be reduced 
to a manageable level by using high-resolution follow-up imaging with adaptive optics 
or the Hubble Space Telescope. Our predictions should help optical transient surveys 
capitalize on their potential for discovering tidal disruption events. 

Key words: galaxies: nuclei — black hole physics 



1 INTRODUCTION 

A massive black hole (BH) at the center of a galaxy 
can tidally disrupt stars that pass within a radius Rt ~ 
-R*(MBH/Af a tar) 1//3 of it, where Mbh and M s t ar are the mass 
of the BH and staiQ, respectively, and R+ is the radius of the 
star. For solar-mass stars, -Rt lies outside the BH's event 
horizon for Mbh <J 1O 8 M0; in these systems, a fraction of 
the stellar debris is expected to flow back towards the BH 
following disruption, releasing a flare of radiation. Although 
the rate of tidal disruption events p er galaxy is uncertain, it 



is likely ~ 10 6 — 10 3 per year (e.g., u . , 
1 19991 : IWang fc Merrittfl2004l ; iDonlev et ai[ f2002) 
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Observing and studying such flares has the potential 
to inform our understanding of accretion physics, the mass 
function of BHs, and the dynamics of stars in the nuclei of 
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galaxies. To date, several candidate events have been discov- 
ered, and many more are likely to be found in the coming 
years. A handful of candidates were discovered in soft X-rays 
by the ROSAT All-Sky Survey and XMM-Newton Slew Sur- 
vey, and several more candidat es have been dis c overed in the 
ultraviolet (UV) by GALEX jKomossal [20021 : iGezari et all 
I20081 . 120091 ); these observations likely probe emission from 
an accretion disk close to the BH. Searches with GALEX 
are ongoing, and several new wide-field, high-cadence opti- 
cal transient surveys have recently started or are planned: 
the Palomar Transient Factory (PTF), Pan-STARRS, and 
the Large Synoptic Survey Telescope (LSST). 

In IStrubbe fc Quataert! l|2009lh we argued that opti- 
cal surveys have the potential to discover tens to hun- 
dreds of tidal disruption events per year. Bright optical 
emission occurs when stellar debris is unbound by the in- 
tense radiation pressure produced by the debris falling back 
to the BH at a super-Eddington rate. We predicted that 
this outflow produces a flare as bright as a supernova, 
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~ 10 43 — few x 10 44 ergs _1 , that lasts for days to weeks. 
These optically luminous flares may be detectable out to 
z ~ 1 (and perhaps beyond). However, optical transient sur- 
veys are finding and will find many bright transients close 
to the centers of galaxies (e.g., active galactic nuclei and 
supernovae) . Identifying tidal disruption events amid an ar- 
ray of far more common transient phenomena thus poses a 
substantial challenge. Observational follow-up is crucial: de- 
tailed multi-wavelength light curves, high-resolution imag- 
ing (to determine that events are truly nuclear), and spec- 
troscopy are all required. In the long term, the last of these 
has the potential to be the most definitive signature of a 
tidal disruption event. 

In this paper, we predict the optical-X-ray spectro- 
scopic signatures of tidal disruption events as a function of 
time, focusing on the outflows produced when the fallback 
rate is super-Eddington. In such outflows, gas outside the 
electron scattering photosphere emits photons and absorbs 
photons from deeper in, producing a spectrum that can con- 
tain emission and absorption features. A separate source of 
spectroscopic features can arise at late times, after the out- 
flow subsides, from the half of the star that gained energy 
upon disruption and is escaping from the BH in the star's 
original orbital plane. The surface of that equatorial ma- 
terial is irradiated by the accretion disk, producing broad 
emission lines (mostly hydrogen) offset in velocity from the 
galactic lines l|Strubbe fc Quataertl 120091 ') . Detecting these 
lines would also be a strong confirmation of a tidal disrup- 
tion event, but they are usually substantially fainter than 
the spectral diagnostics presented here and so will be more 
difficult to observe. 

The remainder of the paper is organized as follows: 
in [JU we review the physics of super-Eddington outflows 
produced during tidal disruption events and describe how 
we use the photoionization code Cloudy to calculate the 
spectral lines of the outflow and how we calculate the line 
profiles. We also critically assess when the assumption of 
thermal equilibrium f or the outflow's emission employed by 
IStrubbe fc Quataertl \200$ ) is valid f i)2.2|l . In Sj3]we describe 
our primary spectroscopic predictions. We then briefly esti- 
mate the rate of supernovae in the nuclei of galaxies (Q, 
since nuclear supernovae are one of the primary sources of 
confusion in optical searches for tidal disruption events. Fi- 
nally, in [J5]we discuss our results and their implications for 
observing and identifying tidal disruption events. 

2 SUPER-EDDINGTON OUTFLOWS 
2.1 Summary of Basic Properties 

We summarize theoretical expectations for the physics 
of tidal disruption events. Following a star's disruption, 
roughly half of the stellar debris becomes bound to the BH, 
falls back to pericenter, a nd shocks; the rate of fallback is 
(|Reedll988l : |Phinnevl[l989l ) 



bit is R p , i? p .3R s = Rp/3Rs (where Rs is the Schwarzschild 



1 M sta 



allback — „ 



3 ^fallback 



V ^fallback / 



-5/3 



where 



(fallback — 20Mg ^ -Rp^Rg 



-3/2, 



(1) 



(2) 



is the period of the most bound debris, the BH mass is 
Mbh = Mq x 1O 6 M0, the pericenter distance of the star's or- 



radius), and the stellar radius 



R^/Re. For M BH <£ 



few x IO 7 M0, the fallback rate predicted by equation |T} 
can be much greater than the Eddington rate M^dd for a 
period of weeks to years; here MEdd = 10LEdd/c 2 , Z/Edd is 
the Eddington luminosity, and 0.1 is the fiducial efficiency 
of converting accretion power to luminosity. 

While the fallback rate is super-Eddington, the stellar 
gas returning to pericenter is so dense that the photons pro- 
duced in the shock are unable to escape and cool the gas; 
in particular, the time for photons to diffuse out of the gas 
is longer than both the inflow time in the disk and the dy- 
namical time characteristic of an outflow. The gas is likely to 
form an adve ctive accretion disk accompanied by powe rful 
outflo ws (e.g.. iKing fc Pounddl20q3l ; lOhsuga et aj]|2005h . 

In IStrubbe fc Quataertl ((2003) , we developed a simple 
model t o describe the outflowing gas (see also related esti- 
mates in lKing fc Pound j2003l ; lRossi fc Begelmanll2009l ). We 
assume the outflowing gas is launched from ~ _Rl = 2R P at 
a rate 



M out = /out M{ allback 

with a terminal velocity 

1>wind = fvV CS c(RL), 



(3) 



(4) 



which is typically 1 — 10% of the speed of light. Radiation 
hydrodynamical simulations of super-Eddington accretion 
show that the density and velocity structure of the outflow- 
ing gas varies with angle, with hig her speed outflows along 
the pole (e.g., lOhsuga et al.l 12005): in the tidal disruption 
context, the outflow properties may also vary with time. To 
account for variations with viewing angle, we consider val- 
ues of /out ranging from 0.01 to 0.3, with / ou t =0.1 as our 
fiducial value, and values of f v ranging from 0.1 to 1, with 
f v = 1 as our fiducial value. We approximate the outflow's 
geometry as spherical, with a density profile 



P(r,t) 



M ou t(t - r/v wind ) 
47rr 2 w win d 



(5) 



inside the outflow where r <J i? e d gc = «wmd(- For r -C R c d gc , 
the density varies as p(r,t) ~ M out (t)/47rr 2 i;wmd. When t ;> 
few x (fallback) the density increases with radius approaching 
R e d S e'- most of the mass is near the edge, within a shell of 
thickness Ar sn cii, so that 



p(Rcdge, t) 



\ /out M s i 



47! -#cdge A '-shcll 



(6) 



Since most of the gas was expelled during a period last- 
ing ~ tfaiiback, Ar s hell is at least ~ Wwindtfallback; to account 
for a possible variation in outflow velocity during the pe- 
riod when most of the wind is launched, we assume that 
the wind speed varies by A« w i n d/«wmd ~ 10%, so that 
Ar sll eU ~ max(u win dtfaiiback, 0.1i; w i n dt). The exact magni- 
tude of Aw w ind is uncertain, but accounting for this sys- 
tematic variation is important because otherwise the shell is 
unphysically narrow and dense at late times. 

At most wavelengths, the dominant opacity in the out- 
flow is electron scattering. The outflow is optically thick to 
electron scattering out to a radius i? pn ,cs, the electron scat- 
tering photosphere, at which R p h,csp(Rph,cs)K cs ~ 1, where 
Kes the opacity to electron scattering: 
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r> ^ , n < f-1 ( fallback \ pl/2 R 

\ M E dd J b 

Because they are trapped by electron scattering, pho- 
tons produced in the shoc k cool adiabatically a s the gas 
expands in the outflow. In IStrubbe fc Quataertl (|2009h we 
assumed that the gas and photons would be in thermal equi- 
librium at the shock so that the outflowing photons would 
have a blackbody spectrum, given by 

vL u = 4K 2 R 2 ph , cs vB„(T ph , cs ) (8) 

where T p h,os is the temperature at the electron scattering 
photosphere. The flow remains supported by radiation pres- 
sure and thus adiabatic expansion causes the temperature 
in the outflow to scale as T a p 1 ^ 3 - In thermal equilib- 
rium, the temperature of the gas and radiation at the shock 
(T eq ) are determined by aT^ ~ tt P rc,gas, i.e., the post-shock 
photon energy density is approximately equal to the bulk 
kinetic energy density of the pre-shock gas. We approxi- 
mate that the gas falls back to pericenter spherically, so 
that Mpre,gas ~ M fa n b ack«Gsc,L / in Rl , where u csc ,L is the es- 
cape velocity of the gas at R^. The temperature at -R p h, es 
in thermal equilibrium is thus 

T ph , cs ~ 2xl0 5 K ( -h-) ( 5/ M- 1/4 ^3C-( 9 ) 

\ /out / \ M E dd / 

After weeks to months, the outflow finally becomes opti- 
cally thin to electron scattering, revealing the accretion disk 
close to the BH. The accretion disk's spectrum is a mul- 
ticolor blackb ody with temperatur e s ~ 1 5 K, described in 
more detail in lStrubbe fc Quataertl (|2009h . Eventually, after 
a time 

t Edd * o.i A'4 2/5 <3R S ™ 3 A*" 3/5 yr (10) 

(where m* = Mstar/M©), the mass fallback rate decreases 
below the Eddington rate, and radiation pressure is no 
longer strong enough to unbind new gas. The previously 
expelled gas continues to expand outwards, becoming a thin 
shell located at ~ R c d S c with a thickness Ar s h e ii and a den- 
sity given in equation (0. From some viewing angles, the 
accretion disk will continue to be seen through this shell, 
with the exact covering fraction of the shell depending on 
the geometry of the outflow at early times, which is some- 
what uncertain. 

2.2 The Applicability of Thermal Equilibrium 

Equations |8j and <(9j assume that the gas and radiation are 
thermally well-coupled from the shock at ~ _R P to the out- 
flow's electron scattering photosphere at i? p h, es - However, 
the post-shock gas and radiation may not have time to reach 
thermal equilibrium before advecting away from the shock. 
Here we quantitatively assess the applicability of thermal 
equilibrium drawing an analogy to the radiation-mediated 
shocks present d uring supernovae and shock-breakout (e.g., 
iKatz et al.|[2O10h . 

In local thermal equilibrium (LTE) , the temperature of 
the gas and radiation are given by T cq , where aT^ ~ u prCig as 
(see the text before eq. [§])■ The dominant continuum emis- 
sion process is free-free emission. The time to reach thermal 
equilibrium £lte is determined by the timescale for free-free 
emission in the post-shock plasma to produce the number 
density of photons required for thermal equilibrium. Prior to 



~ thTE, the gas and radiation are instead in Compton equi- 
librium at a temperature T shock that is substantially larger 
than T eq ijKatz et alJl2D10h . The shock jump conditions im- 
ply that the velocity falls and the density of gas rises by a 
factor of 7 as the gas mo ves from pre - shock to post-shock. 
Using equation (13) from IKatz et alj l|2010l) . we then find 
that the time for the post-shock plasma to thermalize, in 
units of the local dynamical time tdyn = -Rl / iw l , is given 

^ ~ mM ^ s R-X 6 (—) 9/8 (^-) 1B/8 
tdyn 6 p,3Rs \mj \10 day J 

/ D \ 43/16 , j. y 15/ 

~ 0.9M 6 109/4S f #0 m- 9 / 8 r- 27 / 16 f — ({ 

V Rt / V tf allback / 

If £lte < tdyn, the assumption of blackbody emission in 
equations (|8} and ((9JI is reasonable. Equation (|11[) shows 
that this is generally true at early times for Mbh ^ 10 Mq 
and any stellar pericenter distance, and for Mbh <J few x 
10 6 Mq and _R P ~ 3i?s- However, for more massive BHs 
(Mbh ~ 10 7 Mq), the assumption of thermal equilibrium 
is probably poor for t ^ tfaiibackj when most of the stellar 
debris returns to pericenter. Moreover, for a given BH mass 
and stellar pericenter distance, the post-shock plasma is not 
in thermal equilibrium after a time 

+ n k d 47 /30 /?ri*\ 3/5 , 

tnon-thcrm ~ 0.5 Mq' -K p .3 Rs l I day 

Equation (|12|l shows that the assumption of thermal black- 
body emission is likely reasonable for a few weeks for Mbh ~ 
1O 6 M0 and R p ~ Rt- For smaller R p , thermal equilibrium 
breaks down earlier, but events with R p ~ Rt are predicted 
to dominate the rates and a re thus in practice probably the 
most important (Fig. 12 of IStrubbe fc Quataertl l2009h . At 
times t ^ tnon-therm, the temperature at the shock, T a hock, 
is tens to hundreds of keV, scaling as J?"^ for mildly rela- 
tivistic fallback speeds. The radiation emitted at the electron 
scattering photosphere will be cooler than this by a factor 
of ~ (-Rph.es/ Ri,) 2 / 3 due to adiabatic expansion. Compton 
upscattering also likely gives the radiation a power-law spec- 
trum. We will discuss the effects of this non-blackbody emis- 
sion in § H3.1I and 151 but defer a detailed calculation of the 
non-LTE spectrum to future work. 

It is important to note that if the gas and radiation 
are not able to thermalize close to the shock (i.e., for t ;> 
tnon-therm), they are unlikely to thermalize further out in the 
outflow instead: the equilbrium photon number density falls 
with radius, as p, but the equilibrium free-free emissivity 



2 This estimate of £lte i s about 100 times shorter than the most 
naive estimate of the thermalization time, ~ (qjc) -1 , where 
is the Rosseland mean absorption coefficient for free-free interac- 
tions. The time *lte is that associated with waiting for the gas to 
emit photons rather than waiting for it to absorb photons already 
present. The large numerical difference between <lte and ~ a g c 
is primarily due to different weighting of frequencies in the averag- 
ing; in the correct calculatio n, emission of low-frequency photons 
dominates I lKatz et a,l.ll2010f) . while in o^, absorption of photons 
close to the blackbody peak dominates. 
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falls faster, asp 11 / 16 , so there is even less time to come into 
thermal equilibrium at larger radii. 

2.3 Spectroscopic Calculations 

The spectrum of the outflow will be imprinted with spec- 
tral lines produced by the outer layers of gas between the 
photosphere and the edge of the outflow at Rcdgc- The 
photosphere is initially determined by the outflow itself 
(-Rsourcc = -Rph.es) but at later times as the outflow subsides 
and becomes optically thin, the photosphere is set by the 
accretion disk, with -Rsourcc = -Rdisk- The gas outside -Rsource 
absorbs photons released deeper in, through photoionization 
and bound-bound transitions, creating absorption and emis- 
sion features. The outer gas is highly ionized by the central 
source, and maintains photoionization equilibrium so long as 
the recombination time is shorter than the expansion time: 
tree ~ (jieOrcc) -1 < t. Here a Tec is the recombination coef- 
ficient, ~2x 10 _13 cm 3 s _1 for hydrogen at 30, 000 K and 
typically larger for heavier species. 

Early on, most absorption takes place at r ~ R p h, C s, 
where the gas is always in photoionization equilibrium, since 

(-Rph.es) 



[n e (-Rph,es)arcc] 1 
10- 5 ^M 6 5 /X3H S ^^ — 



10 day 



-8/3 

(13) 



At late times, most of the mass in the outflow resides in a 
shell at r ~ -Redge, which dominates the absorption. There, 
the gas falls out of photoionization equilibrium after a few 
months to a few years, since 



tree (-Redge ) 



[n e (R 



edge ) 



3 x 10" 



ft 



3 ^i^R, 



/out 



3/2 



10 day 



(14) 



when Ar s heU ~ Vwindifallbacki t ICC increases even more later 
when Ar she ii ~ 0.1« win dt- 

We determine the ionization and opacity structure of 
the gas outside -Rsourcc by performing photoionization cal- 
culations with version 08. 00 of the publ i cly av ailable code 
Cloudy, last described by iFerland et alj (| 19981 '). To deter- 
mine the spectrum, we then post-process Cloudy's output to 
account for Doppler shifts by the outward motion of the gas. 
We will only show results below for times at which t ICC <J t 
so that photoionization equilibrium is a reasonable approx- 
imation. 

For calculating line profiles, it is useful to divide the 
outflow into two parts: at times t <J tEdd, when the outflow 
is being continuously driven, there is a radially extended 
outflow from -Rsourcc to ~ -Redge- Because -Redge is gener- 
ally much largei0 than -R SO urce, the line-of-sight velocities of 
the gas span a wide range, which causes absorption lines to 
be strongly velocity-broadened. At times t ^ few x tfaiiback, 
there is also a narrow, denser shell at r ~ -Redge, which 



3 In lStrubbe fc Quataerti {2009), we describe an early ~day-long 
phase for small A/bh an d small R p during which -R p h iC s ~ t? c d gc ; 
although there may be interesting spectroscopic features during 
this phase, we focus on later times when the physics of the escap- 
ing photons is more secure. 



contains most of the mass (because most of the mass is un- 
bound at ~ tfaiiback)- In the shell, thermal broadening may 
dominate over velocity broadening. We first qualitatively de- 
scribe the evolution of the spectrum produced by these two 
parts of the outflow, and then explain in more detail how 
we calculate the absorption and emission line profiles. 



2. 3. 1 Three phases of evolution 

For Mbh ^ few 10 7 Mq, the fallback rate is super-Eddington 
and the outflow is optically thick to electron scattering for a 
few weeks to months after disruption. Photons released from 
-Rsource = R-ph,cs with the blackbody spectrum in equation 
© (perhaps with an additional X-ray power-law tail; Q 
photoionize the outer layers of gas between -Rsourcc = -R P h,cs 
and -Redge, which have a density profile given by equation 
(JSJ). This gas spans a wide range in radii, and so produces 
a spectrum of broad absorption lines, whose profiles are de- 
scribed in £12.3.21 when t ;> few x tf a iib ac k, there are also 
narrower absorption lines, described in H2.3.3I The emission 
lines are calculated as described in i|2.3.4l 

At later times, the fallback rate diminishes and the 
electron scattering photosphere of the outflow moves in- 
ward. The entire outflow becomes optically thin to elec- 
tron scattering ~ weeks to months after disruption, but 
it continues to be driven until Mf a iiback falls below MEdd, 
which can be somewhat later (tEdd; eq. I10p . Deep inside 
the outflow, accretion onto the BH proceeds via a thin 
disk, which emits a multicolor blackbody spectrum peakin g 
close to ~ 10 5 K (described in IStrubbe fc Quataerti 120091 , . 
In our calculations for these times, our input spectrum to 
Cloudy is the spectrum of the accretion disk (whose size 
is -Rsourcc = i?disk)- The outflow continues to span a wide 
range in radii, .Rsourcc = -Rdisk <; r <; R c d g e, and so we cal- 
culate absorption line profiles as in c|2.3.2l with additional 
absorption lines from the shell ( H2.3.3P when appropriate. 
The emission lines are again calculated as in c|2.3.4l 

For t > t E dd, Mf a iiback < M E dd and the shocked gas at 
pericenter can cool efficiently; essentially all of the gas thus 
accretes through the disk rather than being blown away. The 
previously expelled material continues to expand out as a 
thin shell, with a density given by equation ([6]) and a radial 
thickness Ar s heli- The shell is irradiated by the blackbody 
emission from the accretion disk. In many cases, the outflow 
is no longer in photoionization equilibrium (see ea. H4[) . but if 
there are places where the outflow velocity is low (/„ ~ 0.1), 
the shell can be dense enough to remain in equilibrium for 
up to several years. In such cases, the absorption lines are all 
narrow f b|2.3.3p . while the emission lines are broad ( >12.3.4|) 
and very faint. 

We now describe how we calculate the line profiles for 
the extended part of the outflow and the shell (see also re- 
lated calculations in ICastor1ll970l . for line profiles in Wolf- 
Rayet star winds). 



2.3.2 Velocity-broadened absorption lines 

Cloudy outputs a table of absorption lines: each entry con- 
tains the line frequency vo, the species (element and ioniza- 
tion stage) producing the transition, and the total optical 
depth Tstat through the (stationary) layer of gas. Cloudy also 
outputs the density distribution ns P ecies('") for each species, 
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and the temperature profile T(r). The stationary optical 
depth is 



Tstat 



species 



V)ao(r) dr 



= (croVth) 



^-species 



Source ^h^) 
ccics i 



dr 



(15) 

(16) 
(17) 



where ao(r) is the cross section of the transition at line cen- 
ter and Vth(r) is the thermal velocity of the gas. Because 
the lines provided by Cloudy are thermally broadened, the 
quantity (aoVth) is independent of radius. 

For simplicity, we assume that the gas flows out radially 
with a spatially and temporally constant velocity, v — u w j n d 
(eq. SJ, superposed by small thermal motions. Since the gas 
is optically thin at most frequencies, its temperature reg- 
ulates to T ~ 10 4 — 10 5 K, leading to thermal velocities 
Wth ~ 10 — 30kms _1 <C n w i n d ~ 0.01c — 0.1c. Because the 
photoionizing source — the electron scattering photosphere 
or accretion disk — is spatially extended, its radiation orig- 
inates from impact parameters 6 ranging from (center of 
the source) to -Rsource (edge of the source). At a given impact 
parameter b, our line of sight passes through gas moving at 
projected line-of-sight velocities «los = w-windyl — {b/r) 2 , 
where r ranges from -Rsource to -Redgc. 

A transition of frequency vq can absorb photons of rest 
frequency v at places in the wind where ulos satisfies v = 
1/0(1 + vlos/c). These locations are centered at radii 

r abs = b (18) 

V 1 - ("LOs/tAvind) 2 

with a small spread along the line of sight, Al, due to ran- 
dom thermal motion of the gas: 



Al = r abs 



Vth(r a bs) \ fr, 



f wind 



(19) 



Thus, for a given rest frequency and impact parameter, the 
optical depth to a given transition ifl 



T b ~ (n s peciesO"oA^)| r 



^stat^abs ^species 



(20) 



& 2 Vwind/s 

Approximating the source as a spherical isotropic emitter, 
we find the overall transmitted power: 



rt: 



r/2 



2e~ Tb sin 9 cos 9 d9 , 



(21) 



where sin 9 = b/R SOUICO . When multiple transitions con- 
tribute absorption at frequency v, we replace Tb above with 
the sum of their corresponding optical depths. We account 
for continuous absorption processes as well; these are gen- 
erally important only for hv ;> 0.3 keV since the gas is so 
highly ionized. 

2.3.3 Thermally-broadened absorption lines 

At times t J> few x tfaiiback, the bulk of the previously ex- 
pelled gas forms a thin shell at r ~ -Redgc, as described above 

4 We use Tstat and -/species because Cloudy does not output line 
opacities as a function of radius. 



in H2.ll While velocity broadening is the dominant broad- 
ening mechanism for absorption at -Rsource <J r <^ -Redge, the 
line-of-sight velocity's variation with impact parameter at 
r ~ -Redge may become less than t> t h after at most a few 
months if the wind velocity is close to constant in time. For 
example, at times when -Rsource = -R P h,os, 



A«los 
«th 



0.03 



R 



ph.es 



Re, 



I «wind \ 



/out , f 16/3 „ll/2 
~~ff M S H P,3Rs 



~rj\ 



30 kms 



10 day 



-16/3 



.(22) 



(The variation in «los along the line of sight through 
the narrow shell is even smaller than this, by a factor of 
Ar s heii /-Redge-) Consequently, the random thermal motion 
may dominate the broadening of absorption lines produced 
in the shell; such lines are therefore narrow, with a linewidth 
Av = vo(vth/c), and are blueshifted to ~ fo(l + t>wind/ c ). 

If the wind speed varies in time by more than ~ 0.1%, 
velocity broadening instead dominates over thermal broad- 
ening, leading to wider and shallower absorption lines. The 
wind speed variation Aw w i n d could be substantial, e.g., ~ 
O.lfwind, but is highly uncertain, so we consider thermal 
broadening in the shell as a lower limit. 

We run separate Cloudy calculations for the narrow 
shells, which give the thermally-broadened optical depth 
Tth(^) as a function of which we blueshift by w w ind/c. 
At times t <J iEdd, we multiply the velocity- broadened spec- 
trum described above by e -7 **'"'; at later times, we multiply 
the continuum produced at -Rdisk by this factor. 



2.3.4 Emission lines 

The gas at r ^ -Rsource can also produce emission via ra- 
diative recombination and radiative decay of (collisionally 
or radiatively) excited atoms/ions. Because of the large ve- 
locities in the outflow, the gas is effectively optically thin 
to all of these photons, even at the energies of resonance 
lines. As a result, we can observe emission from gas hav- 
ing line-of-sight velocities ulos ~ — Vwind up to «los ~ 
+«wind(l — R^sourcc/R^cdgc) 1 ' 72 ; we approximate this with a 
Gaussian line profile centered on the rest energy of the line, 
and slightly truncated on the red side. This approximation 
is based on analogy to Monte Carlo calc ulations of emis- 
sion from expanding Lyman alpha blobs |Verhamme et all 
2006). We treat emission lines produced in the narrow shell 
the same way, since line photons emitted from essentially 
any part of the shell can reach us. 

The large velocities in the outflow imply that the re- 
sulting emission lines are extremely broad; in most of our 
calculations, the emission is so spread out that the lines 
will be orders of magnitude fainter than the transmitted 
spectrum. The overall emission vL^ s then consists of faint 
continuum and very broad emission lines centered on the 
rest frequencies of the transitions. The observable spectrum 
is the sum of this emission and the transmitted light from 
equation (|21|l . i.e., 

(23) 



Lout t trans . r emis 



Although the physical processes are similar in tidal 
disruption events and broad absorption line quasars (BAL 
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QSOs), we note that tidal disruption events have substan- 
tially smaller emission line equivalent widths; this compari- 
son is discussed further in ij5] 

3 PREDICTED SPECTRA 

We now use the methodology of 321 to calculate spectra as a 
function of time due to the disruption of a solar-type star, 
varying the BH mass and pericenter distance of the stellar 
orbit. We assume solar abundances. We focus on solar type 
stars because they are among the most abundant stars at 
~ 1 — 10 pc in galactic bulges, which is where most of the 
disrupted stars originate. The three fiducial models we con- 
sider are: M BH = 1O 6 M and R p = 3R S ; Mbh = 1O 6 M 
and R p = .Rt; and Mbh = 1O 7 M and R p = Rt- Our fidu- 
cial model for the outflow takes f v — 1 and / out =0.1, but 
later we vary these values. 

To start, we assume that the outflow produces a thermal 
blackbody spectrum that photoionizes the surrounding gas; 
in i|3.1l we consider the effects of (harder) non-thermal emis- 
sion on the predicted spectra. The thermalization time esti- 
mate in >I2. 2 1 implies that thermal equilibrium in the outflow 
is maintained for t <; 0.5 and 10 days for Mbh = 10 6 M , 
R p = 3Rs and Rp = Rt, while it fails at later times; the 
outflow is never in thermal equilibrium for A/bh = 1O 7 M0. 
Given, however, the uncertainties in the precise thermaliza- 
tion time we show thermal outflow models for a range of 
timescales. 

Figure [1] shows our predicted spectra at various times 
after disruption, from the optical to the extreme ultraviolet 
(EUV). Significant lines in the far ultraviolet (FUV) to opti- 
cal are labeled in the top and bottom panels. The top panel 
of Figure [2] shows a zoomed-in version of the 1000 - 2000 A 
region of the spectrum for our fiducial Mbh ~ lO 6 A/ and 
i?p = Rt model. For the latter, we plot the luminosity den- 
sity L\ normalized to its maximum value over the wave- 
length interval 1000 — 2000 A and vertically offset different 
curves for clarity. 

In Figures[T]and[2] we show results from t ~ tfailback, the 
peak of fallback and outflow, to when the outflow shell falls 
out of photoionization equilibrium^ (The outflow's pho- 
tosphere -Rph,es — where most absorption takes place while 
the fallback rate is super-Eddington — is always in photoion- 
ization equilibrium; see beginning of ^2.31 ) Photoionization 
equilibrium in the shell typically fails when, or somewhat 
before, the fallback rate reaches the Eddington rate. More 
precisely, it fails at ~ 30 days for A/bh = 10 6 Mq, R p = 3i?s, 
/out = 0.1; 100 days for Mbh = 1O 6 M , R p = Rt, 
/out = 0.1, and 200 days for Mbh = 1O 6 M , R p = R T , 
/out = 0.3; while fEdd ~ 30, 400 and 400 days for these 
three models respectively. For Mbh = 10 M©, R p = -Rt, 
the shell of gas at ~ R c d S c is never in photoionization equi- 
librium, and so we do not include any contribution from the 
shell. 

The spectral features visible in Figure [T] are exclusively 
absorption lines; the emission lines are so broadened by the 
large range in line-of-sight velocity that they become un- 
detectable above the blackbody continuum. The outflow as 

5 For Mbh = 1O 6 M , R p = 3Rs (top panel), the earliest time 
we depict is t = 3 days because -R pn ,os ~ Hedge until ~ 1 day (see 
footnote [3j . 



a whole and the dense thin shell at -R c d gc typically both 
contribute to the absorption lines. For A/bh = 1O 6 M , 
R p = 3-Rs, the shell dominates the absorption lines at 
A ^> 800 A; the shell and outflow as a whole both contribute 
at a wide range of wavelengths for 1O 6 M , R p — Rt; and 
the shell does not contribute for 1O 7 M , R p = Rt because 
it is never in photoionization equilibrium. The absorption 
lines from the extended part of the outflow are strongly 
blueshifted and typically (while t < tEdd) very broaqj, 
with linewidths of ~ lAvmd ~ 0.01 — 0.1 c. Superimposed 
on these broad lines are the thermally-broadened narrow 
lines produced in the outer shell of gas, with linewidths of 
~ 10- 3Qkms _1 . 

Most of the absorption lines are in the UV (A £ 2000 A; 
hv ;> 10 eV), with few features in the optical. The reason for 
this has two parts. First, consider the ionization parameter 
Uhv for a species whose ionization energy is £i on = hv: be- 
cause the gas density is relatively low while the incident 
spectrum is luminous and peaks at energy fei/ poa k ~ 10 eV 
or higher, t/ioev 1 and so species having £i on <^ 10 eV are 
almost fully ionized. (For example, the ionization parame- 
ter for hydrogen is typically (7h ~ 10 3 — 10 5 , and hydro- 
gen's neutral fraction is typically 10 -8 — 10~ 10 .) Secondly, 
most atoms/ions are in the ground state: the radiative decay 
rates are fast compared to the rates of photoionization and 
collisional excitation. As a result, almost all of the species 
present in the flow have £i on ^ hvp ca ,k <; 10 eV, and mem- 
bers of those species are in the ground state. Transitions 
from the ground state have energies similar to £i on , which 
is J> 10 eV, so most spectral lines have energies ^ 10 eV, 
i.e., in the UV rather than optical. In particular, some com- 
mon prominent lines in the UV at A > 1000 A are He II 
(1640 A), C IV (1548 + 1551A), Si IV (1394 + 1403A), O 
IV (1400 + 1401 + 1405 + 1407A), N V (1239 + 1243A), Ly- 
man a (1216A), and O VI (1032 + 1038A). These are similar 
to the lines observed in the spectra of BAL QSOs — this is 
not surprising given that the physical conditions are simi- 
lar. When the density is high enough for optical lines to be 
present (because U and hv pca .k are lower), these optical lines 
are mostly the lines of H I and He II. 

Figure [T] shows that the minimum energy /maximum 
wavelength of the absorption lines in the optical-UV de- 
pends both on time and on the parameters of the tidal dis- 
ruption (e.g., AfsH and R P ). For example, the typical ab- 
sorption lines shift to shorter wavelengths at later times. The 
physical origin of these dependencies can be understood as 
follows. For all times and outflow parameters that we con- 
sider, Uhv in the outflow is 2> 1 for hv <J /t^ poa k, and so 
the approximate minimum energjQ of significant absorption 
lines present in the spectra is ~ /i^ pea k/few. So long as the 
outflow is optically thick (most times depicted), the peak 
energy of the incident spectrum is set by the temperature 
of the electron scattering photosphere T phies . We can thus 



6 The sawtooth shape of the absorption lines at early times in 
FigurefTJis a result of our assumption of spatially constant velocity 
(see i|2,3.2l l; a more realistic velocity gradient would shift some 
absorption from ~ ^^(l + fwind/c) towards hvo, leading to a less 
abrupt change in the spectrum at ~ hv<^{\ + ^> w ind/c). 

7 High-n lines of hydrogen at wavelengths of tens of microns are 
also typically optically thick due to /-mixing collisions with ions. 
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Figure 1. Predicted spectra for our three fiducial tidal disruption flares at several different times after disruption. Labeled tickmarks 
identify the blueshifted positions of the strongest long-wavelength lines. A zoomed in view of the FUV region of the middle panel is in 
Figure [2] At early times, the continuum emission in these calculations is produced by a super-Eddington outflow while at later times it 
is produced by the accretion disk close to the black hole. These calculations assume that the outflow is able to thermalize completely 
f i|2.2t ; Fig. [4] shows results for incomplete thermalization. 



understand the variation in the typical energy of spectral 
lines by considering the scalings for T p h, C s in equation ([9]). 
For example, at a fixed time after disruption, a model with 
Mbh = 1O 6 M0, R p — 3Rs has a hotter electron scattering 
photosphere than a model with Mbh = 1O 6 M0, R p = Rt- 
This is why the spectra for the former situation have fewer 
low-energy lines in Figure[TJ Mbh = 10 7 Mq, R p = Rt is yet 
hotter at fixed time, and so has even fewer longer wavelength 
lines. 



These arguments also help to explain the time evolution 
of the spectra in Figure [T] During the optically thick phase 
(^ I2.3.1|) . the continuum radiation produced by the outflow 
becomes harder with time as R p h,as moves inward, while the 
luminosity remains high. As a result, the lines present in the 
spectrum tend to have shorter wavelengths (higher energies) 
at later times. Once the outflow becomes optically thin, the 
continuum spectrum is ~ 30 eV — 100 eV emission from the 
accretion disk, and most absorption lines have hv ^ 30 eV 
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Figure 2. Predicted spectra varying the mass-loss rate in the outflow (via / ou t ; eq. [3J , focusing on the wavelength region 1000 — 2000 A. 
Each spectrum is portrayed as L\ , normalized by its maximum value on that wavelength range, with an added integer constant to offset 
curves for clarity. Labeled tickmarks identify the blueshifted positions of the lines. Top panel: Mbh = 10 6 A/q, R p = Rt, / ou t = 0.1, 
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stronger for larger mass outflow rates (larger /out), because the continuum emission has a lower temperature (T p h, e s). 



and remain very broad ( i|2.3.2[) . In some cases, narrow lines 
produced in the shell ( i]2.3.3[) are superimposed on these 
broad lines. The lines would become purely narrow after 
tEdd because nearly all of the outflow is then in a thin shell, 
but for our fiducial outflow parameters, the outflow falls out 
of photoionization equilibrium before t^dd- 

As mentioned in H2.ll numerical simulations suggest 
that the velocity and density of the outflow will vary with 
latitude, with higher spe ed outflows along the pole rela- 
tive to the equator (e.g., lOhsuga et alj [20051 ). To consider 
how the spectrum of a tidal disruption event may vary with 
viewing angle, Figure [3] shows spectra for different values 
of /out and f v (eqs. [3] and [J]). Significant NUV/optical 
lines are labeled; the lower panel of Figure [2] highlights the 
A = 1000 - 2000 A part of the spectrum for the / ou t = 0.3, 
f v = l model. Because the density at r ~ -Redge is larger if 
the outflow is slower (smaller /„) the shell remains in pho- 
toionization equilibrium longer for the models with /„ =0.1 
in Figured for 4000 days (Mbh = 10 6 M Q , R p = R T ) and 
1000 days (M BH = 1O 7 M , R p = R T ). 

The ionization parameter Uhv is large for all of these 
variations about our fiducial models. It is thus again the tem- 
perature of the continuum radiation T p h, os that determines 
the approximate minimum energy of the spectral lines. For 
/out = 0.3, Tph.cs is lower than for / ou t = 0.1, and so there 
are more and deeper FUV-optical lines when / ou t is larger 
(compare the top panel of Figure with the middle panel 
of Figure [T] or the two panels of Figure [2]); for / ou t = 0.01 



(not shown), there are virtually no lines with A ;> 1500 A. 
Similarly, T p h, C s is lower for lower outflow velocities. For 
Mbh = 1O 6 M0, R p — Rt, fv = 0.1, there are in fact many 
optical absorption lines (mostly hydrogen Balmer and He II) 
early on; for Mbh = 10 7 Mq, R p — Rt, fv = 0.1, there are 
many lines at ~ 1000 — 2000 A, though no lines in the op- 
tical. Additionally, Ha and He II A6560 (blended together) 
can be seen in emission at late times for Mbh = 1O 8 M0, 
i?p = -Rt, fv = 0.1. Since the outflow velocity is lower, the 
density at r ~ R c d gc is larger and so recombination happens 
more frequently, leading to emission lines; furthermore, the 
lower outflow velocity produces less broadening and so the 
emission lines are brighter above the blackbody continuum. 

3.1 Implications of an X-ray Power-law 

The candidate tidal disruption events discovered in the 
ROSAT All-Sky Survey and GALEX Deep Imaging Survey 
all show so ft X-ray spectra jKomossal |2002| ; iGezari et al.1 
|200S| . |2009| ). However, by analogy to the observed spectra 
of AGN, it is possible that some tidal disruption spectra 
will contain a high-energy power-law tail extending from the 
peak in the blackbody continuum to hard X-rays. Such an 
X-ray power-law component might also be produced as a re- 
sult of incomplete thermalization at the shock at pericenter 
where matter falls back to the BH ( N2.2[l . To consider the 
observational effects of such X-rays, we carried out Cloudy 
calculations using an input spectrum that consists of the 
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Figure 3. Predicted spectra showing the effects of varying the mass-loss rate (via /out; eq[3]l and outflow speed (via f v \ eq. [4} in the 
super-Eddington wind. Labeled tickmarks identify the blueshifted positions of the strongest long-wavelength lines. Slower and/or higher 
mass-loss rate winds have softer continuum emission and are thus more likely to produce optical or near-UV lines (compare these results 
with Fig.0. Top panel: A/ BH = 1O 6 M , R p = R T , / out = 0.3, f v = 1, t = 30 days, 100 days; middle panel: M BH = 10 6 M Q , R p = R T , 
f out = 0.1, /.„ = 0.1, t = 30 days, 100 days, 300 days, 1000 days; bottom panel: M BH = 10 7 M Q , R p = R T , f out = 0.1, /„ = 0.1, t = 30 
days, 100 days, 300 days, 1000 days. A zoomed in view of the FUV region of the middle panel is in Figure [2] 



blackbody spectrum described in i]2.3.1l plus a power-law 
tail that has 10% of the blackbody luminosity, a photon in- 
dex of 3 {vL y oc u" 1 ), and that begins at the frequency 
where its emission equals the blackbody emission. (Such a 
spectrum may be appropriate for partial thermalization; in 
a moment, we consider the even harder spectrum expected 
from Compton equilibrium if the shock is not thermalized.) 

Figure [4] shows the 10 keV to 1 /im spectra for several of 



our models including this X-ray power-law: Mbh = 10 Mq, 
R p = R T , fv = 1 (top panel); M BH = 10 7 M o , R p = R T , 
f v = 1 (middle panel); and M BH = 10 6 , R p = R T , fv = 0.1 
(bottom panel). Due to the different peak energies of the 
blackbody components, the actual power-law luminosity 
above 1 keV varies for the different models as a function 
of time, with uL v (> lkeV)/L bo i = 10" 4 - 2 x 10~ 2 (top 
panel), 6 x 10~ 3 - 10~ 2 (middle panel), and 5 x 10~ 5 - 10~ 2 
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Figure 4. Predicted spectra, including the presence of an X-ray power-law tail with photon index T = 3 that carries 10% of the 
blackbody luminosity (appropriate when the shock at pericenter does not have time to thermalize completely). Optical and UV lines are 
mostly absent because the gas is so highly ionized (cf. pure blackbody continuum in Figs, [l] and although there is typically significant 
absorption in the EUV and soft X-rays. When T = 2, there are even fewer lines and little continuum absorption. 



(bottom panel). These models are thus reasonably conserva- 
tive in terms of the contribution of the X-ray emission to the 
bolometric luminosity. Nonetheless, the presence of the hard 
X-ray emission significantly changes the resulting spectra. 

The hard incident spectrum photoionizes the gas to a 
higher degree than the pure blackbody incident spectrum. 
Species that have ionization energies £j on < 10 eV are thus 
more scarce, and so the optical absorption lines seen in the 



previous section disappear (except for Mbh = 10 Mq, R p = 
Rt, f v = 0.1, where Ha still has optical depth ~ 1). Most 
of the FUV lines disappear as well; the Lya line and a few 
others can still be faintly visible early on for Mbh = 1O 6 M0, 
especially if f v ~ 0.1. Fe XXI A1354 is sometimes the longest 
wavelength UV/optical line. 

So long as the outflow is optically thick to electron scat- 
tering, the X-ray tail can show many absorption features. 
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These include both continuum absorption and individual 
absorption lines. Figure [4] shows that for Mbh = 10 6 Mq, 
R P = Rt, there is a deep continuum absorption trough, ex- 
tending from ~ 1 keV up to ~ 5 keV; for f v =0.1 the outflow 
is denser and the trough extends down to ~ 0.3 keV. This 
feature is somewhat weaker for Mbh = 1O 6 M0, R p = 37?s 
and absent for Mbh = 1O 7 M0. There are also many absorp- 
tion lines superposed on the power-law tail and absorption 
trough, provided by highly ionized Ar, Ca, Fe, Mg, Mn, Ni, 
and Si (among others). The specific lines and line strengths 
vary significantly between the different models, so it is dif- 
ficult to predict exactly which lines will be in the spectrum. 
As the electron scattering photosphere moves inward and 
becomes hotter, Figure [4] shows that the X-ray luminosity 
actually increases in time given our assumption of a fixed 
power-law with a photon index of 3. Once the outflow be- 
comes optically thin to electron scattering and the accretion 
disk provides the incident spectrum with hi/pc^ ~ 0.1 keV 
(typically after a few months), the gas is so highly ionized 
that most of the absorption and emission features in the 
spectrum disappear. 

If instead of vL v oc v~ , the X-ray power-law tail has 
a flatter spectrum with a photon index of 2 (extending up 
to ~ 100 keV), the gas is even more highly ionized. There 
are thus fewer X-ray lines (though still a significant number) 
and little continuum absorption. There are no lines at all in 
the FUV/optical (apart from very faint Lya and Fe XXI 
A1354 for /„ = 0.1). 



4 SUPERNOVA RATES IN GALACTIC 
NUCLEI 

Tidal disruption flares are expected to be similar to super- 
novae in their overall brightness and timescale. Supernovae 
thus represent a significant source of contamination when 
trying to discover and study tidal disruption flares in op- 
tical transient surveys. In particular, tidal disruption flares 
may be confused with supernovae that appear coincident 
with the galactic nucleus within the spatial resolution of the 
observations. To quantify this source of confusion, we esti- 
mate the rate of Type II and Type la supernovae that take 
place within a distance R res of a galaxy's nucleus. To do so, 
we first estimate the disk and bulge stellar mass within R TCS 
as a function of BH mass. 

For a BH of m ass Mbh, the stell ar mass of the bulge is 
M* lB ~ 700M BH l|Haring fe Rbj2004h . We estimate the stel- 
lar mass of the disk using the ratio of bulg e (B) to total mass 
(T) from Figure 3 of iHopkins et all (|2009h . which is based on 
data from lBalcells et alj \200H ): if M, R < 10 10 M Q , B/T ~ 
0.1; for larger masses, B/T ~ 3 log(M», s /lO lo M ) - 2.9 
(where total mass M»,t = disk mass M»,d + bulge mass 
M^b). We determine the half-li ght radius of the d isk Ri/2,d 
using the results from SDSS in lShen et alj l|2003l ) (their eq. 
18) and convert to scale radius via 7? SC aic,i3 = -Ri/2,d/1-68. 
We then find the mass of the stellar disk inside R ICS by as- 
suming that the surface density profile is exponential: 



-Rrcs/-R SC alc,D 

M*,zj(< R res ) ~ M,, D I xe~ x dx. 



(24) 
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Figure 5. Supernova rates within 0.5" (black) and 0.05" (blue) 
of the galactic nucleus at z = 0.1, for Type II supernovae (solid) 
and Type la supernovae (dashed); 0.5" is the typical astrometric 
accuracy of ground-based transient surveys (J. Bloom, personal 
communication) while 0.05" is the spatial resolution of HST or 
ground-based adaptiv e optics. If the rate of tidal disruption per 
galaxy is 10 — 5 yr" _1 jPonlev et alj l200l red dot-dashed line), 
then nuclear supernovae will outnumber tidal disruption events 
by about two orders of magnitude for ground-based optical tran- 
sient surveys. For HST/AO, the nuclear supernova rate will be 
comparable to the tidal disruption rate, reducing the contamina- 
tion dramatically. 



Rm.B ho st ing a given BH usingj equation (17) in 
Sh en et al.l (|2003l ) and convert it to the scale radius of 
the de Vaucouleurs profile using -R SC aic,B = i?i/ 2 ,s/3461 
l|Binnev fe MerrifieldlTl99ct ). We then find the mass of the 



bulge inside R TCB using 



M,, B (< Rr, 



M, 



20160 



■Rres /-R S ca 



dx . 



(25) 



Similarly, we calculate the half-light radius of the bulge 



The rate of Type II supernovae depends on the instan- 
taneous star formation rate (SFR) in the galaxy. For galaxies 
that are actively forming stars, the galaxy-integrated SFR 
as a function of galaxy stellar mass (M* = M*,xj + M*,b) is 
log SFR = 0.671og(M« /10 10 M w ) - 6.19, where SFR is in 
M Q yr _1 l|Noeske et al.l |2007v FI. Assuming that the spatial 
distribution of star formation tracks stellar mass, the SFR 
within R rcs is SFR(< R [cs ) ~ [M*, D (< i? rcs )/M,, D ] x SFR. 
The nuclear Type II supernova rate is then 



8 Note the erratum in IShen et"aD l|2003h 's Table 1, so that b = 
2.88 x 10~ 6 . 

9 We extend the relation down to M* ~ 1(PMq, and we sup- 
pose that the relation is proportional to (1 + z) 3,4 , using the 
re dshift dependence of the volume-integrated star formation rate 
in lYiiksel et alj 1120081 ) and in agree ment with the approximate 
redshift dependence reported bv lNoeske et al. 
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r n (< R T , 



w - 2fD SFR«R 



M yr- 



res I —1 



(26) 



To account for the fact that only disk galaxies typically form 
stars, equation (|26[1 includes a multiplicative factor given by 
the fraction of disk galaxies at the given M», 



fo 



driE /din M» 
dntot /d In M* 



(27) 



where dn/d In M* is the num ber density of e l liptica l galaxies 
(E) or all galaxie s (tot) from lBernardi et al.l (|2010l ') localhPI 
and iDrorv et ail l|2009l . their Table 3) at higher redshift. 
With this factor of fo, our estimate of the supernova rate 
per galaxy statistically takes into account that some systems 
of a given M* are already passive and thus will not have 
significant numbers of Type II supernovae. 

Following recent galaxy integrated results, we estimate 
the rate of nuclear Type la supernovae given both the nu- 
clear stellar mass and the nuclear star formation rate: 



r Ia (< Rr, 



A 



[M„, B +M*, D ){< Rr, 



+ B[fn 



1O 1O M 

SFR(< R les ) 
lOMgyr- 1 



(28) 



yr 



where A = 4.4 x 10 yr~ 
|Scannapieco fc BildsterJl2005l 1. 

The typical accuracy to which ground-based optical 
transient surveys can determine the location of a transient 
is ~ 0.5" (J. Bloom, personal communication), which corre- 
sponds to R ICS = 0.9 kpc at z = 0.1 and R rcs = 4.1 kpc at 
a = l. Figure [5] shows the Type II and Type la supernova 
rates per galaxy within 0.9 kpc of the galactic nucleus, as 
a function of BH mass; this choice of _R rcs corresponds to a 
ground-based survey observing at z ~ 0.1. Predicted tidal 
disruption rates depend on BH mass and galaxy structure 



and typically range from 



10" 



10 3 yr 1 per galaxy. 



Using candid a te de tections in the ROSAT All-Sky Survey, 
iDonlev et all (|2002h estimated a rate of ~ 10 J tidal dis- 
ruptions per year per galaxy. For the latter, Figure [S] shows 
that nuclear supernova rates are typically several orders of 
magnitude larger than tidal disruption rates. High resolu- 
tion photometry with the Hubble Space Telescope (HST) or 
ground-based adaptive optics (AO) can decrease the nuclear 
SN rates to few x 10 -5 yr per galaxy, reducing confusion to 
order unity. We discuss the implications of these estimates 
of nuclear supernova rates in the next section. 

5 DISCUSSION 

We have calculated the spectroscopic signatures of outflows 
produced by super-Eddington accretion during the tidal dis- 
ruption of stars by massive black holes. Although there 
are some uncertainties in the continuum emission, mass 
outflow rate and kinematics, we find a number of reason- 
ably robust conclusions: 1) the spectrum will show strong 
absorption lines that are blueshifted relative to the host 
galaxy, 2) the absorption lines will typically be very broad 



10 Paramet ers for mass functions of ellipticals and total are from 
Table B2 in Bcrnardi ct al. defining ellipticals by a con- 

centration index larger than 2.86 as recommended by these au- 
thors. 



(~ 0.01 — 0.1c), though often with a narrow thermally- 
broadened (~ 30kms _1 ) component that dominates as the 
outflow subsides and becomes a thin shell at later times, 3) 
if the continuum spectrum is largely a blackbody, the ab- 
sorption lines will be most prominent at UV wavelengths 
(e.g., C IV, Si IV, O IV, N V, Lyman a and O VI). In addi- 
tion, if there is a lower velocity component to the outflow, 
which is plausible b ased on simulations of radiatively ineffi- 
cient accretion fe.g-. lQhsuga et aLll2005T ). we find that there 
will also be optical absorption (and possibly emission) lines, 
in particular Ha, H/3, He II A6560, and He II A4860 (Fig. [3]). 
The optical/UV lines will, however, be largely absent from 
the spectrum if the continuum emission is harder or contains 
an X-ray power-law tail (either because of non-thermal pro- 
cesses or incomplete thermalization) — the gas would then 
be too highly ionized. In that case, the dominant absorp- 
tion lines are in the soft X-rays (Fig. [4]). This highlights the 
importance of X-ray observations of tidal disruption events 
coeval with optical and UV spectroscopic observations, in 
order to properly interpret the presence or absence of opti- 
cal/UV lines. 

Having summarized our key results, we now describe 
several uncertainties in our spectroscopic predictions. Our 
calculations assume a blackbody spectrum of photons re- 
leased from the electron scattering photosphere of the out- 
flow, with an 'optional' phenomenologically motivated X-ray 
power-law tail. We find, however, that when gas returning 
to pericenter shocks at ~ 2R P , the thermalization between 
the photons and gas is a strong function of BH mass, stellar 
pericenter distance, and time since disruption f i|2.2[) . The 
thermalization is often incomplete. This would harden the 
continuum spectrum, eliminating many of the optical-UV 
spectroscopic signatures of tidal disruption events. More- 
over, the prominent optical continuum emission predicted by 
IStrubbe fc Quataerj (|2009h from tidal disruption outflows 
requires reasonable thermalization in the post-shock plasma. 
Such thermalization is the most likely for tidal disruptions 
around lower mass BHs with Mbh <J few 1O 6 M0 and at early 
times after disruption, t <J 2 weeks (eg. Ill[) . This estimate 
assumes spherical fallback and is thus somewhat conserva- 
tive, since stellar debris is focused into a thin stream as it 
falls back after disruption (|Kochaneklll994r i ; a higher density 
at the shock would lead to a longer period of thermal black- 
body emission, perhaps up to a month for Mbh ~ 1O 6 M0. 
Nonetheless, these results emphasize the importance of high 
cadence observations with rapid follow-up in optical searches 
for tidal disruption flares. In the future, observing the spec- 
trum of a tidal disruption event change as the outflow falls 
out of thermal equilibrium would provide strong constraints 
on the physics of radiation-dominated shocks, which are im- 
portant in other astrophysical environments such as shock 
break-out in supernovae. 

Another uncertainty in our spectroscopic predictions is 
related to the mass loading and speed of the outflow (pa- 
rameterized by / ou t [eq. [3] and f v [eq. U in our calcula- 
tions), and the outflow's geometry. Radiation hydrodynamic 
simulations of super-Eddington bla ck hole feeding sugge st 
reasonable values for / out and /„ : lOhsuga et al.l (|2005T ) 's 
simulation at lOOM Edd obtains f o u t ~ 0.1 and /„ ~ 1 
(our fiducial values), lohsugal (|2007l Vs similar simulations of 
feeding at M Edd - 300M Bdd imply / out ~ 0.1 - 0.4 (de- 
pending on viscous parameter a) and outflow velocities of 
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~ 0.1e-0.3e; rTakeuchi et all (|2009h 's similar simulations at 
100M E dd and 300M Ed d imply / out - 0.8. iDotan fc Shavivl 
calculate a super- Eddington accretion model in which 
instabilities make the gas inhomogeneous, and find that 
/out ~ 0. 5 - 0.7 for feed i ng ra tes 5A'/ E dd - 20M B dd- We 
note that lLodato fc Rossi l|201ll ) calculate tidal disruption 
light curves us ing our super-Ed d ingto n outflow model (eqs. 
[7] - [9]) with IDotan fc Shavivl (|2010l ) 's /out results and a 



more detailed model for Mfaiiback(t) ( Lodato et al. 2009) 



they find similar results to those in IStrubbe fc Quataertl 
(2009) (optical luminosities are similar, and emission from 
the outflow lasts a factor of few times longer) . In we used 
fiducial values of / ou t ~ 0.1 and /„ ~ 1 and showed the re- 
sults of varying those parameters: larger mass outflow frac- 
tions lead to more absorption lines at longer wavelengths, 
and slower outflow velocities lead to narrower absorption 
lines, the presence of optical lines, and even the presence of 
some emission lines. In reality, the mass-loading and out- 
flow kinematics may vary with viewing angle; the geometry 
of the outflow is uncertain, and we d o not capture this ef- 
fect w ith o ur simple spher i cal m odels. lOhsuga fc Mineshigel 
i|2007l) and lTakeuchi et all (|2009l ) find outflow opening (half) 
angles of ~ 30 degrees, suggesting that from some viewing 
angles, accretion disk and outflow may both be visible. 

5.1 Observational Prospects 

The prominent UV lines predicted here are challenging to 
observe for several reasons. First, the UV emission can be 
obscured by dust, particularly along lines of sight through 
the host galaxy's disk; however, observations of similar ab- 
sorption lines in BAL QSO spectra suggest that at least 
some lines of sight will have low obscuration. Second, obser- 
vations at extreme UV wavelengths (1000 A <J A 100 A), 
where most of the predicted spectral features lie, must take 
place from space and are technically difficult. Encouragingly, 
the Space Telescope Imaging Spectrograph and the Cosmic 
Origins Spectrograph aboard the Hubble Space Telescope 
(HST) should be able to observe spectroscopic features like 
those predicted here in the far UV. 

The optical to X-ray spectra predicted here apply to 
tidal disruption flares having early-time super-Eddington 
outflows, which requires Mbh <J few x 1O 7 M0. To detect 
such flares in the first place using optical transient surveys, 
significant contaminants such as variable active galactic nu- 
clei (AGN) and supernovae must be excluded. Luminous 
AGN fueled by other means are ;> 10 3 times more common 
than tidal disruption events and could in principle produce 
optically bright 'flares.' Typically, however, AGN show opti- 
cal emission lines in their spectra (like the Balmer lines and 
[O III]). We have shown that during the super-Eddington 
outflow phase, tidal disruption events are unlikely to show 
such optical lines; furthermore, if tidal disruption events do 
show optical lines, Balmer and He II lines are the only rea- 
sonable candidates (perhaps also very faint He I and S IX), 
because the densities are too high for collisionally excited 
lines like [O III] which almost always appear in AGN spec- 
tra. These conclusions about Balmer versus [O III] lines also 
hold duri ng the later phase of a tida l disruption event (de- 
scribed in lStrubbe fc Quataerj[2009l ) in which the accretion 
disk irradiates the unbound stellar debris. 

BAL QSOs (a subset of quasars) are more physically 
and spectroscopically similar to tidal disruption events, both 



involving a bright central continuum source driving an out- 
flow, which gives rise to blueshifted UV absorption lines. 
However, BAL QSOs show strong emission lines while tidal 
disruption events likely will not, for th e follo wing two-part 
reason, which draws on iMurrav et all jl995h 's theoretical 
work on BAL QSOs. First, the tidal disruption outflow is 
typically much more ionized than a BAL QSO wind, so op- 
tical depths for resonance lines in tidal disruption events 
are at most comparable to, and are often much less than, 
optical depths in BAL QSOs. (The optical depth to true 
absorption — rather than scattering — by resonance transi- 
tions is typically much less than one in both cases; the emis- 
sion is thus effectively thin and so the smaller optical depth 
for tidal disruption outflows implies less emission.) Secondly, 
BAL QSO winds are thought to originate well outside the 
source of continuum radiation (the gas is not as highly ion- 
ized at larger radii, so UV resonance lines are optically thick 
enough for radiation pressure to drive the wind). Thus, al- 
though we see the hot continuum source through the cooler 
wind, the wind can produce significant emission lines (in 
addition to strong absorption) because the wind's emitting 
area is much larger than the area of the continuum source 
jMurrav et al.lll995h . In tidal disruption events, by contrast, 
those emitting areas are typically the same: the region where 
most line emission occurs and the outer edge of the con- 
tinuum source are both typically at the electron scattering 
photosphere. This, along with the small true optical depth, 
implies that tidal disruption outflows are expected to show 
little or no line emission. The situation may be more com- 
plicated if the outflow does not fully thermalize, as there 
could be additional line emission inside the electron scat- 
tering photosphere that we have not considered here(3 In 
addition to this important difference in the presence of emis- 
sion lines, note that BAL QSOs are generally associated with 
higher-mass BHs (Mbh <; 1O 8 M0), which cannot tidally dis- 
rupt solar- type stars outside the horizon. 

Another significant contaminant in optical searches for 
tidal disruption events is supernovae that take place in the 
nuclei of galaxies, because the luminosities and timescales 
are similar, and because supernovae are initially blue like 
tidal disruption events. Spectroscopy will help: supernovae 
typically show a forest of deep optical absorption lines, quite 
unlike our predictions for tidal disruption events. Type Iln 
supernovae may still pose a particular challenge: the super- 
nova plows into surrounding circumstellar medium, shock- 
heating the gas to high enough temperatures to suppress 
the optical lines. However, these supernovae begin to cool 
after several months, becoming redder and producing opti- 
cal emission lines at late times. Tidal disruption events, by 
contrast, become hotter with time and are not expected to 
produce similar optical emission lines. 

To quantify how much nuclear supernovae will 'contam- 
inate' searches for tidal disruption events, we estimated the 
rate of supernovae within the spatial resolution (~ 0.5") of 
ground-based optical transient surveys (^4] and Fig- [5j : at 
z ~ 0.1 the nuclear Type II rate is a few x 10~ 3 yr -1 per 
galax y, of which ~ 10% are probably Type Iln (|Li et al.l 
2010); the Type la rate is a factor of 2 — 3 smaller than 



11 If the outflow runs into relatively dense circumnuclear gas, 
additional line emission may be possible there as well. 
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the Type II rate. These rates of Type II and Type Iln su- 
pernovae are ^ 2 and 1 order(s) of magnitude greater than 
the tid al disruption rate per galaxy inferred bv lDonlev et al.l 
(2002) using ROSAT. The ROSAT constraints, however, are 
largely on more massive black holes (which are both more 
luminous and more likely to be prominent in the X-rays) 
so it is possible that the tidal disruption rate is different in 
the lower mass systems most likely to produce optically lu- 
minous emission from super-Eddington outflows. The high 
rate of nuclear supernovae estimated in Figure [5] emphasizes 
the importance of high-resolution follow-up imaging with 
HST or adaptive optics, which can reduce the rate of nu- 
clear supernovae 'false positives' in optical tidal disruption 
searches by a factor of ^ 50 (Fig. [5]l . 

Although the super-Eddington outflows produced dur- 
ing the tidal disruptio n of stars are in principle r eadily de- 
tectable out to z ^ 1 (jStrubbe &: Quat aert 2009), restrict- 
ing candidate events to hosts with z ~ 0.1 is likely a good 
strategy for minimizing interference from supernovae: the 
rate of nuclear supernovae increases rapidly with redshift, 
to ~ 10 _1 yr _1 per galaxy at z ~ 1. This increase is due 
to both the increase in star formation at high redshift and 
the increasing fraction of a galaxy that lies within the point- 
spread function of the observation. 

The fiducial outflow model in IStrubbe fc Quataert! 

(2009) (/out = 0.1; f v — 1; disruption rate per galaxy = 
10 _5 yr _1 ) predicts that ~ 8 tidal disruptions per year can 
be detected at z ^ 0.1 for a survey like the Palomar Tran- 
sient Factory. Relatively nearby events have the additional 
advantages that they are less expensive to follow up spectro- 
scopically and it is easier to characterize their host galaxies. 
For these comparatively nearby events, we predict that the 
detection probability in optical surveys is relatively indepen- 
dent of Mbh for Mbh few x 10 7 M Q , so it should help to 
restrict candidates to bulges less massive than ~ 1O 1O M 
(the weak dependence of the detection probability on Mbh 
is also important because tidal disruption flares from lower 
mass BHs ~ lO 6 A/ are the most likely to be in thermal 
equilibrium and thus to have prominent optical continuum 
emission in the first place). Restricting follow-up to systems 
with little ongoing star formation and/or old stellar popula- 
tions would help further minimize the number of nuclear su- 
pernovae. Observationally, the star formation rate in galax- 
ies, and hence the supernova rate, is relatively bimodal, and 
thus some galaxies will have nuclear supernova rates (par- 
ticularly type II rates) smaller than estimated in Figure [5] 
This is especially true for more massive BHs, but ~ 10% of 
lower massive systems (M bh <^ 1Q 7 M ) are li kely to be rela- 
tively passive as well fe.g.. lBernardi et al.ll201oT ). If these can 
be identified in advance (via, e.g., prior Sloan Digital Sky 
Survey observations), they may well be the most promising 
systems in which to follow up nuclear transients; this selec- 
tion would, however, decrease the predicted detection rate 
to ~ 1 yr~ at z <J 0.1 for a survey like PTF. 

Given the recen t plethora of luminous supernovae (e.g., 
iQuimbv et al.l l2009). it is unclear how well single-band pho- 
tometry alone will be able to distinguish tidal disruption 
events from nuclear supernovae. If this proves difficult, as 
we suspect is likely, the spectroscopic predictions for tidal 
disruption e vents presented here , and the color evolution 
predicted in IStrubbe fc Quataert! l|2009T l (the outflow pho- 
tosphere becomes hotter with time — becoming bluer if the 



observing band is close to the peak, or showing no color evo- 
lution if the observing band is on the Rayleigh- Jeans tail), 
may prove particularly useful for identifying and character- 
izing tidal disruption flares. 

5.2 Optically-selected candidates 

The first two optically-selected tidal disruption candidates 
were announced shortly after the submission of this pa- 
per, found in Stripe 8 2 of the Sloan Digital Sky Survey 
(van Velzen et all 120101 1. These candidates have observed 
properties broadely consistent with our predictions, with 
observed g-band luminosities of ~ 10 43 ergs _1 and no color 
evolution. The optical data for candidate "TDE1" is reason- 
ably approximated by our model for the super-Eddington 
outflow due to the disruption of a solar-type star by a BH 
of mass ~ 1O 7 M0 at R p ~ _Rt, with / ou t ~ 0.1 [eq. [3] and 
f v ~ 0.1 [eq. H]. In detail, the optical data for candidate 
"TDE2" is harder to approximate with our simple model: 
although the luminosity can be reproduced by a relatively 
large / ou t or relatively small /„, the optical colors and only 
gentle fading of the event are less consistent with our model. 
These properties may be mo re consistent with a model like 
that of lLoeb fc Ulmerl (|l997f ) in which falling-back gas set- 
tles into a steady hydrostatic atmosphere rather than be- 
coming unbound in a true outflow; this may be appropriate 
if the fallback rate is never highly super-Eddington. The 
large optical luminosity of these two candidate events sug- 
gests that outflows may have relatively low velocities, and is 
encouraging for the optical detection of future tidal disrup- 
tion candidates. 



ACKNOWLEDGMENTS 

We thank Josh Bloom, Kevin Bundy, Phil Chang, Wei- 
dong Li, Norm Murray, Kristen Shapiro, and Nathan Smith 
for helpful conversations; we also especially thank David 
Strubbe for help with C++ calculations. Support for EQ was 
provided in part by the Miller Institute for Basic Research 
in Science, University of California Berkeley and the David 
and Lucile Packard Foundation. Support for program num- 
ber AR 12151 was provided by NASA through a grant from 
the Space Telescope Science Institute, which is operated by 
the Association of Universities for Research in Astronomy, 
Inc., under NASA contract NAS5-26555. 



REFERENCES 

Balcells M., Graham A. W., Peletier R. F., 2007, ApJ, 665, 
1084 

Bernardi M., Shankar F., Hyde J. B., Mei S., Marulli F., 

Sheth R. K., 2010, MNRAS, 404, 2087 
Binney J., Merrifield M., 1998, Galactic astronomy, Binney, 

J. & Merrifield, M., ed. 
Castor J. I., 1970, MNRAS, 149, 111 

Donley J. L., Brandt W. N., Eracleous M., Boiler T., 2002, 

AJ, 124, 1308 
Dotan C, Shaviv N. J., 2010, ArXiv e-prints 
Drory N., Bundy K., Leauthaud A., Scoville N., Capak P., 

Ilbert O., Kartaltepe J. S., Kneib J. P., McCracken H. J., 



© 2011 RAS, MNRAS OOO.rnfl5l 



Spectroscopic Signatures of the Tidal Disruption of Stars 15 



Salvato M., Sanders D. B., Thompson D., Willott C. J., 

2009, ApJ, 707, 1595 
Ferland G. J., Korista K. T., Verner D. A., et al., 1998, 

PASP, 110, 761 
Gezari S., Basa S., Martin D. C, et al., 2008, ApJ, 676, 

944 

Gezari S., Heckman T., Cenko S. B., et al., 2009, ApJ, 698, 
1367 

Haring N., Rix H., 2004, ApJL, 604, L89 

Hopkins P. F., Somerville R. S., Cox T. J., Hernquist L., 

Jogee S., Keres D., Ma C, Robertson B., Stewart K., 2009, 

MNRAS, 397, 802 
Katz B., Budnik R., Waxman E., 2010, ApJ, 716, 781 
King A. R., Pounds K. A., 2003, MNRAS, 345, 657 
Kochanek C. S., 1994, ApJ, 422, 508 

Komossa S., 2002, in Rev. Mod. Ast., Vol. 15, , pp. 27-+ 
Li W., Leaman J., Chornock R., Filippenko A. V., Poznan- 
ski D., Ganeshalingam M., Wang X., Modjaz M., Jha S., 
Foley R. J., Smith N., 2010, ArXiv e-prints 
Lodato G., King A. R., Pringle J. E., 2009, MNRAS, 392, 
332 

Lodato G., Rossi E. M., 2011, MNRAS, 410, 359 
Loeb A., Ulmer A., 1997, ApJ, 489, 573 
Magorrian J., Tremaine S., 1999, MNRAS, 309, 447 
Murray N., Chiang J., Grossman S. A., Voit G. M., 1995, 
ApJ, 451, 498 

Noeske K. G., Weiner B. J., Faber S. M., Papovich C, 
Koo D. C, Somerville R. S., Bundy K., Conselice C. J., 
Newman J. A., Schiminovich D., Le Floc'h E., Coil A. L., 
Rieke G. H., Lotz J. M., Primack J. R., Barmby P., Cooper 
M. C, Davis M., Ellis R. S., Fazio G. G., Guhathakurta 
P., Huang J., Kassin S. A., Martin D. C, Phillips A. C, 
Rich R. M., Small T. A., Willmer C. N. A., Wilson G., 

2007, ApJL, 660, L43 
Ohsuga K., 2007, ApJ, 659, 205 

Ohsuga K., Mineshige S., 2007, ApJ, 670, 1283 
Ohsuga K., Mori M., Nakamoto T., Mineshige S., 2005, 
ApJ, 628, 368 

Phinney E. S., 1989, in IAU Symposium, Vol. 136, The 
Center of the Galaxy, M. Morris, ed., pp. 543 — h 

Quimby R. M., Kulkarni S. R., Kasliwal M. M., Gal- Yam 
A., Arcavi I., Sullivan M., Nugent P., Thomas R., How- 
ell D. A., Bildsten L., Bloom J. S., Theissen C, Law N., 
Dekany R., Rahmer G., Hale D., Smith R., Ofek E. O., 
Zolkower J., Velur V., Walters R., Henning J., Bui K., 
McKenna D., Poznanski D., Cenko S. B., Levitan D., 
2009, ArXiv e-prints 

Rees M. J., 1988, Nature, 333, 523 

Rossi E. M., Begelman M. C, 2009, MNRAS, 392, 1451 
Scannapieco E., Bildsten L., 2005, ApJL, 629, L85 
Shen S., Mo H. J., White S. D. M., Blanton M. R., Kauff- 
mann G., Voges W., Brinkmann J., Csabai I., 2003, MN- 
RAS, 343, 978 
Strubbe L. E., Quataert E., 2009, MNRAS, 400, 2070 
Takeuchi S., Mineshige S., Ohsuga K., 2009, PASJ, 61, 783 
van Velzen S., Farrar G. R., Gezari S., et al., 2010, ArXiv 
e-prints 

Verhamme A., Schaerer D., Maselli A., 2006, A&A, 460, 
397 

Wang J., Merritt D., 2004, ApJ, 600, 149 

Yiiksel H., Kistler M. D., Beacom J. F., Hopkins A. M., 

2008, ApJL, 683, L5 

© 2011 RAS, MNRAS 000,rrHT5l 



